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ABSTRACT 

The problems associated with data acquisition and reduction in the U.Va. 
superconducting magnetic suspension and balance facility are similar to 
those in free-flight ranges (or tunnels). The model undergoes a "Quasi-six- 
degree-of- freedom" motion whjch must be monitored both in position and angu- 
lar orientation from which the aerodynamics must be inferred. The data 
acquisition problem is made more difficult because geometric constraints 
prevent direct visual access to the model in the Mach 3 wind tunnel. The 
methods, accuracies and problems associated with the acquisition of data 
are discussed. 
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1.0 Introduction 


The problems associated with data reduction in the U.Va. superconduct- 
ing magnetic suspension and balance facility (SMSB) are similar to those in 
free-f light ranges (or tunnels). A complete description of this facility 
can be found in reference I. The model undergoes a "quasi -si x-degree-of 
freedom" motion which must be monitored in both position and angular orien- 
tation, from which the aerodynamics must be inferred. The advantages of 
this fac i I i ty over a convent iona I free f I ight fac i I ity is the ability to 
test over "long" times and gather more data, making the determination of 
model aerodynamics more accurate. In addition the increased length of time 
enables the transient part of the motion to decay allowing observation and 
analysis of steady-state motion. This promises to be a useful technique, 
as wi I I be described below. 

The precision with which we can determine the aerodynamics depends on 
two factors: first, the accuracy of the data acquisition technique used; 

and, second, the manner in which the errors propagate through the mathe- 
matical motion model and data reduction techniques. Here we will primarily 
be concerned with the latter; however, since it is an integral part of the 
data reduction scheme we will first discuss the data acquisition problem. 

2.0 Data Acquisition 

There are three methods for data acquisition in the U.Va. magnetic sus- 
pension facility. These are shown in Figures I, 2, and 3. The first of 
these-the optical sensor (Figure I), used to provide feedback for the con- 
trol system, will also provide the primary data source for the data reduc- 
tion process. The optical system is a conventional light beam-photocell 
system designed and calibrated to give position and angular data. The model 
geometry is one of the optical sensor design criteria, and hence, changing 
model geometry may require a new optical sensor design and calibration. 

The optical system will hopefully be replaced by an electromagnetic 
sensor of the MIT type (Figure 2) which at this time is still under develop- 
ment. The main problems with the electromagnetic system arise from its use 
near a high alternating current source. The electromagnetic sensor will re- 
quire less area in the annulus around the tunnel and thus allow a larger 
diameter test section. The primary motivation for development of the 
electromagnetic sensor is its invariance to model changes. 2 

The fiber optic system^ (Figure 3) was designed primarily for use as a 
visual cue for the operator, who, due to tfre helium dewar, does not have a 
direct line of sight to the model. The distortion created (mostly barrel 
type), shown in Figure 4 can be compensated for by an extensive calibration 
procedure. This source of data is difficult to i ncorporate into the data 
reduction process since it requires a relatively large preprocess i ng effort 
and its estimated accuracy is an order of magnitude less than the other two 
systems as is illustrated in Table I. 


Table I 

Data Acquisition Systems Accuracies (estimates) 


Sensor Type 


Position 


Ang I es 


Optica I 

Electromagnetic 
Fiber Optic 

3.0 Control Technique 


. I mm (3 axes) 
. I mm (3 axes) 
I mm (3 axes) 


.05 degrees (2 planes) 
.01 degrees (2 planes) 
. I degree (I plane) 


In order to understand the concepts of data reduction as applied to 
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the model motion in the U.Va. magnetic suspension facility it is first 
necessary to understand the "quasi -six-degree-of-f reedom" nature of the 
motion. By "quas i -six-degree-of-f reedom" we mean the model is free to both 
rotate and translate at frequencies above some cutoff imposed by the SMSB 
control' system. The feedback controller is designed to control only low 
frequency (0(10 hz)) and DC components of the model motion leaving the high 
frequency motion "untouched". 


The maximum excursion of the model 
frequencies given by 


o 


max 


m ax 
mi. 


from the tunnel centerline at high 


< 3 . 0 - 1 ) 


where a^x is the maximum angle of attack, l z the pitch moment of inertia, 
m, the mass, and l the distance between center of pressure and center of 
gravity. 

For a typical model with a natural aerodynamic frequency of 35 hz the 
maximum displacement can be kept within I cm of the centerline. 

In principle a controller which will leave the model aerodynamics "un- 
touched" (i.e. not affect the roots of the characteristic equation asso- 
ciated with the aerodynamics) can be designed. It is one that requires 
feedback in position, velocity, angle and angular rate, the gains of each 
being determined by the method of Bass and Gura.4 This method requires the 
knowledge of the aerodynamic properties of the model a priori, and freedom 
to use feedback in all the problem variables. Although feasible, this 
method is less desirable than a simple position-velocity controller. 

An analysis of a simple position-velocity controller with a 3 hz natural 
frequency and /~2 damping ratio has been carried out based on rather crude 
estimates of model aerodynamics. The indications are that there is little, 
if any, interaction with the model motion due to aerodynamics. 

Two sources of error were examined: 

1. Errors due to uncertainties in model aerodynamics. 

2. .Errors due to uncertainties in the position of the magnetic center 
with respect to the center of mass. 

As can be seen from Figure 5 there is an insignificant effect on the 
damped natural frequency of the model and a small effect (about 2%) on the 
damping exponent. The error introduced by the controller into the damping 
exponent can easily be calibrated with wind off and used to compensate the 
later results of the inversion process to obtain aerodynamics. For the 
steady-state case a frequency response analysis for a typical model with 
and without controller was conducted. The results are given in Table II. 

Table I I 


Frequency Response Analysis 


Driving Frequency 
(Rad/Sec) 


(Rad/Ft- 1 b ) 


l X/M ex+l 

(Cm/Ft-lb) 


Free 

Control led 

Free 

Control 1 ed 

150 

1 . 56036 

1.56981 

17.4107 

17.5269 

175 

1 .91601 

1 .92581 

15.7085 

15.7964 

200 

2.45552 

2.46566 

15.4143 

15.4838 

225 

2.96322 

2.96845 

14.698 

14.7285 

250 

2.65051 

2.64787 

10.6494 

10.6414 

275 

1 .88953 

1 .88687 

6.27438 

6.26686 

300 

1.33748 

1 .33608 

3.73196 

3.7287 
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The conclusions +o be drawn from the effects of the controller on the 
characteristic equation and the frequency response of the model are as 
f o 1 1 ows : 

1. For transient analyses the damping exponent may be affect by some 
small percentage (about 2-5?). This can be compensated for a 
posteriori . 

2. There is insignificant effect on the frequency. 

3. The forced steady-state motion is essentially unaffected by the 
control ler. 

Thus for the analysis of aerodynamics discussed in the next sections 
the model will be considered to be in free-f light with no inputs due to the 
feedback control system. This, it is felt, is the unique feature of this 
wind tunnel system - "long term" free-f light data. 

4.0 Mathematical Models 

4.0.1 The first mathematical model to be considered is well known , 
linearized equations for a rolling missile with trigonal or greater symmetry 


[2yD - D - ip rt (C 7 + > D)]A-C2y'+ C„ -ip^C, UDA = F gx+ (4.0. 1-1) 


Z Z. r OZ„ Z i z r oz 

a a pg pg q pr 


-[C + C D + ip (C D)]A+Ci 0 D 2 - C -ip (i. - C )]DA = M , (4.0. 1-2) 

u m m- K o m i B m K o A m - 1 ext 

a a pg q pr 


where the C's are the aerodynamic coefficients, y the nondimensional mass, 
i^ and ig nondimensional inertias, p Q the nondimensional roll rate, D the 
derivative operator, A the complex angle of attack, A the complex orienta- 
tion angle, F ex -|. and M ex+ external driving functions. These equations have 
the familiar quadricyclic solution for either of the variables A or A, e.g. 


( A i + iu>i)t (X 2 + i<i»3t 

A = Kje + K 2 e + K 3 e + K 4 


(4.0. 1-3) 


where Kj is the initial amplitude of the mode, Xj is the damping rate of the 
mode, coj is the frequency of the mode, and is the phase angle of the mode. 
The subscripts I, 2, 3, and 4 refer to the precession, nutation, rolling 
trim, and nonrolling trim modes respectively. 


The constants Kj , X;, o>f , and <j>j contain the information needed to 
obtain the aerodynamic coefficients. The precision to which the coefficients 
can be determined depends on the precision to which the data is known; 
examples of this will be given below. 


Equation 4.0. 1-3 contains the information needed to fit both transient 
(all four modes of motion) and steady-state (just the K3 and K4 modes of 
motion) data. The inversion process for transient case yields the aero- 
dynamics in one run using the following relationships 

(4.0. 1-4) 


C = (oj 1 o) 2 - XjX 2 )2 I J pU 2 Sd 

a 


c + c = C (x 

mq m* 


y 

pUS 


1 + * 2) + 2m 


( " C z " VU4I /pUSd 2 
a y 


m 


P3 


wiX 2 + io 2 X^ 

^ + (J ) 2 


pUS 

2m 


(-C - C n )34l /pUSd 2 

z. D x K 
a 


(4.0. 1-5) 


(4.0. 1-6) 


where the drag coefficient and I ift curve slope must be obtained by other 
means. The drag coefficient, Cq, determination is straightforward - being 
proportional to the force required to prevent the model from moving along 
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the tunnel axis. This force is easily determined to a high level of accuracy 
from the currents in the coils. The slope of the lift coefficient, C 2 (as 
well as other translational aerodynamics) can be obtained from a a 

standard swerve reduction program5 (proportional to the lateral distance 
traveled ) . 

The useful part of the application of the quadricycl ic solution to both 
transient and steady-state data lies in the ability to write either the real 
or imaginary part of the solution; e.g. 

3 - Re(A) = K 1 cos(a)it + <}>!) + K 2 cos(u) 2 t + $ 2 ) + l^cosc^t + RetKt*) (4.0. 1-7) 

and still have all the information contained. This enables the application 
of the techniques indicated to data from a single plane. 

4.0.2 The second model to be considered is the fitting of the obser- 
vations to the equations of motion using a technique which we will call the 
"Brute-Force" method. 7 For an axi symmetric model only the z force and M 
pitching moment equations are necessary for the inversion; however, data on 
all kinematic variables is needed: 


wC + qC t pvC t DwC + pDvC + prC = y(Dw - qu + pv) 
Z w Z q z pv Z Dw Z pDv Z pr 


(4. 0.2-1 ) 


wC + qC + pvC + DwC + pDvC 
m ^ m r m m n r 

w q pv Dw 


m t prGm = igDq - <l c - i A >pr 


pDv 


pr 


+ i F (p 2 - r 2 )-i c (qr - Dp) + i n (qp - Dr) (4. 0.2-2) 


where v, w, p, q, r are the nond imensiona I kinematic velocities and the i ? s 
are moments and products of inertia. This model is fitted by reduction to 
a set of algebraic equations as described below. 


4.0. 3 The third model is a specialization of the first one (eq. 4.0. 1-2) 
Here only steady-state motion is considered, perhaps the most unique model 
for^free flight facility to be using. The advantage to steady-state reduc- 
tion is, of course, the increased accuracy to which the data can be deter- 
mined, having many cycles of data to "smooth" over. 

4.0. 4 The last model considered is the full six degree of freedom 
equations of motion given in reference 5. Here the motion is allowed to 
include nonlinear aerodynamics as well as nonlinear inertia terms. 


5.0 Data Reduction Techniques 

Three conceptually different classes of methods have been investigated 
for extracting aerodynamics from available observations. These three classes 
are referred to here as 

1. Differential Correction Methods, 

2. "Brute Force" Method, and 

3. Steady State Analysis Method. 

Most conventional procedures belong to class (I). Our particular adaptations 
are in some aspects unique, as will be explained below, but in general 
represent the state-of-the-art of this approach. Classes (2) and (3) are 
original approaches g$sewing out of our research efforts at U.Va. 

5.0.1 Differential Correction Methods 


Here we are referring to the class of numerical methods which success- 
ively improve preliminary values for the unknown parameters in a given 
mathematical model until the computed output agrees with observations in 
some optimum sense (in our case, minimizing the weighted sum of squares of 
observed-minus-computed residuals). We have employed two differential 
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( 5 . 0 . 1 - 1 ) 


correction formulations In our data reduction analyses, these are 
AP = (A T WA)“ 1 A T W At, 
and 


where 


AP = -(~.) 1/2 G 
G G 

AP = n x I matrix of corrections to parameters, 

Af i m x I matrix of observation residuals, 

A = m x n Jacobian matrix of partial derivatives of the 
m observables with respect to the n parameters, 
evaluated with the current parameter estimates, 

W = m x m weighting matrix, 

<j> = AIRWAY, 

® T - 1 * * • H- 1 

3 Pi c a P n c 


( 5 . 0 . 1 - 2 ) 

( 5 . 0 . 1 - 3 ) 
( 5 . 0 . 1 - 4 ) 

( 5 . 0 . 1 - 5 ) 
( 5 . 0 . 1 - 6 ) 
( 5 . 0 . 1 - 7 ) 

( 5 . 0 . 1 - 8 ) 


and T 

AC = AP AP, AC assigned empirically. (5.0. 1-9) 

The reader is referred to reference 8 for theoretical derivations and 
discussions of (5.0. 1 — I ) and (5.0. 1-2). 


The first differential correction formula (5.0. 1-1) is the classic 
least squares solution. The second formula (5.0. 1-2) is the method of 
gradients ("steepest descent") solution for minimizing an arbitrary function. 


Evaluation of the derivative matrix (5.0. 1-5) is often a source of 
numerical difficulty. For analytic algebraic observation equations, we 
have developed, and used extensively a computer program which completely 
automates the process of partial differentiation. This process was employed 
with the quadricycl ic solution as given in 4.0.1. For those cases in which 
the full six-degree-of-freedom equations (4.0.4) were integrated, we adopted 
a process known as parametric differentiation for computation of the elemenls 
of the observation Jacobian (5.0. 1-5). This procedure8 develops a set of 
m x n differential equations (one for each element of A) which can be 
integrated simultaneously with the equations of motion. These equations 
follow from straight forward partial differentiation of the equations of 
motion. 


Comparing the method of gradients correction equation (5.0. 1-2) with 
the least square correction equation (5-0.1— I ), we note that use of (5.0. 1-2) 
eliminates the necessity of inverting (A^WA), but introduces the necessity 
of controlling convergence rate by logically assigning AC (5.0. 1-9) in 
(5.0. 1-2). Our experience indicates that (5.0. 1-2) is a valid alternative 
to (5.0. I— I ), but should be employed only in the event that (A^WA) is so 
poorly conditioned that numerical inversion is impossible. We have found 
that the classical least square solution is typically an order of magnitude 
more efficient as the basis for least square differential corrections. 

5.0.2 The "Brute Force" Method 

Equations 4.0. 1 — 1 , 2 are solvable by differentiating the data numeri- 
cally (using a five point central differencing scheme) and assuming all the 
observable kinematics to be known quantities. Thus the equations of motion 
are reduced to a set of linear algebraic equations in the aerodynamics which 
can be inverted to obtain the aerodynamics. As one might suspect, the 
accuracy of this method is highly sensitive to errors in observed data, since 
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numerical differentiation is being performed. This method also requires data 
on both position and angles in two orthogonal planes, however it is capable 
of handling a more sophisticated Ci.e. nonlinear) model of the aerodynamic 
forces and moments. 

For more details on this "Brute Force" method see reference 7. Compari- 
son of this method with the others will be presented below. 

5.0.3 Steady-State Analysis 

5.0.3. I Equations 

As stated before, the coupled, comp I ex, second order, linear differen- 
tial equations which describe the motion of our model reduce to two complex 
algebraic equations when only the steady-state response is considered. 

The resulting equations can be arranged^ in the following fashion: 

2y iuiA 


Rea 1 

^ A ~ 


o 

Imag 

2y iu>A 
o 

Real 



o 

Imag 

M ' 

‘tt’ + : 

o 


) = Cp (C 
K 0 2 


+ C 


)Ju) 


"pr 


+ C 


PS 

)3w 


+ Cc 2 2 

a 


+ C-p c 

K o z 


P6 


B- 


o A- 


m 


2y 


(5.0.3. 1-1 ) 


(5.0.3. 1-2) 


(5.0.3. 1-3) 


(5.0.3. 1-4) 


P6 


q a 

The terms appearing on the left hand side of the equations are all observa- 
bles (or in the case of the last equation - computable before they are 
needed). These quantities are determined for each of several frequencies 
on the frequency response curve. Due to the periodic non-damped nature of 
the steady state solution A 0 and A 0 can be determined using a simple least 
square procedure or fourier analysis to obtain amp I itude, phases, and 
frequencies. Since a linear model has been assumed, observations of one 
plane of data (both angle and velocity) is sufficient. Each of the first 
four equations are valid for the n-points used on the frequency response 
curves. Therefore, we have n sets of equations whose solution is a simple 
non-iterative least squares reduction for the coefficients. 

The major advantages of this method are its simplicity (no interation 
necessary) and its relative insensitivity to reasonable measurement errors. 

The determinable coefficients include all of the coefficients on the 
right hand side of equations 5.0. 3.1-1, 2, 3 and 4. This, it should be 
noted, i nc I udes inertia terms. 

5. 0.3. 2 Magnetic Investigation of Resonance 

In the steady state case the use of an oblate spheroid for the 
support element will also allow for an investigation of a resonance curve. In 
the previous section this was shown to be sufficient to determine the models 
aerodynamics. Use of an oblate spheroid provides an additional "spring constant" 
term (which is proportional to the magnitude of the main field) in the rota- 
tional equation of motion. 

For discussion's here we assume no translation, the models motion is 
given by: 
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{-[C +C + C D + ip (C D)3 + i n D 2 - C D - 1 p ( i a D — C D) > = Me 
mm. mi r o mi B m . r o A rn _ o 

PP cj P" 


iuit 


“ a bal “ 


C5.0.3.2-I) 


iwt 


Assuming the steady state solution (A = A,-,e and for simplicity assume 
P Q = 0, then rearranging and separating the equation into real and imaginary 
parts; 

M 

Real (-£) = -C + C-i R a ) 2 -Cl 
& m . d m 

° a bal- 

M 

Imag (-j2.) = (C + C )u> 

3 A mm- 

o q a 


(5. 0.3. 2-2) 


For numerous runs, all with the same u> but different C 


m 


(5. 0.3. 2-3) 

the problem 


a BAL 

becomes the same as that described in the previous section. It should be 


noted that the "variable" is now C 


m. 


not "to" as before. 


*bal 

This method will be investigated numerical ly in the near future, no 
results are as yet available. 


A disadvantage of this method that is that the translation equation of 
motion is unchanged. All the quantities appearing are constant, making it 
necessary to vary to (as well as C ) to evaluate all of the aerodynamic 
coefficients involved. "^BAL 

6.0 Numerical Results 

A comparative numerical study of how observational statistics propagate 
through the data reduction methods into statistics of the determined aero- 
dynamic derivatives has been carried out. Observations were simulated by 
corrupting perfect (computed) values of the observables by adding Gaussian 
random relative errors. Several noise samples were taken at each noise 
level (a); and each of the several applicable data reduction techniques 
were employed to determine the correspond i ng values for the aerodynamic 
derivatives. From these results, small sample statistics of the determined 
derivatives were computed for each method. Typical results of these analyses 
are displayed in figures 6, 7, 8 and 9 for a 15° included angle cone. 

All the data are presented as percent standard deviation of the aero- 
dynamic coefficients versus percent random, noised superimposed on the data. 
The major observation to be made is the consistent superiority of the steady- 
state method over the others. The translational derivatives (not shown) 
follow the same pattern with the steady-state method yielding the most 
accurate inversion at a given noise level. The errors noted in the 
differential correction methods are approximately the same as other investi- 
gators using these methods have found them to be. 

One interesting fact to be reported on in detail in a future publication 
is the ability to separate C m and C m> for reasonable noise levels using the 
"Brute Force" method. ^ ® 


t Curves I and 2 have a slightly higher positional noise level than the 
others, however our experience indicates that the moment coefficients 
are not extremely sensitive to positional noise. 
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7.0 Conclusions 


An analysis of several methods for obtaining aerodynamic coefficients 
from the U.Va. superconducting magnetic suspension and balance wind tunnel 
system has been carried out. The method for inverting steady state free- 
f light motion yields more precise aerodynamic coefficients than transient 
methods at the same measurement noise level. 
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Figure 4 Distorted Cone 



Figure 5 Root Locus of Model Aerodynamics 




